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Abstract: The influence on the I-V characteristics and tunnel magnetoresistance 
(TMR), of impurities embedded into the insulating barrier (I)separating the two 
ferromagnetic electrodes (F) of a magnetic tunnel junction, was theoretically inves- 
tigated. When the energy of the electron's bound state at the impurity site is close 
to the Fermi energy, it is shown that the current and TMR are strongly enhanced 
in the vicinity of the impurity. If the position of the impurity inside the barrier is 
asymmetric, e.g. closer to one of the interfaces F/I, the I-V characteristic exhibits a 
quasidiode behavior. The case of a single impurity and of a random distribution of 
impurities within a plane were both studied. 

Magnetic tunnel junctions (MTJ), consist of two metallic ferromagnetic electrodes sepa- 
rated by an insulating barrier. They typically exhibit tunnel magnetoresistance (TMR) of 



the order of 50% associated with a change in the relative orientation of t 



le magnetization in 
especially due to 



the two ferromagnetic electrodes. They attract a lot of attention 
their applications in several spin-electronic devices especially in non-volatile MRAM (Mag- 
netic Random Access Memory). In a pioneer paper 0|, a theory of TMR for ideal MTJ 
(without defect) was presence of 

different types of defects within the barrier can dramatically affect the I-V characteristics 
and TMR amplitude. In these papers, the current, averaged over the cross-section of the 
system, was calculated. However, it is also interesting to investigate the local current den- 
sity and TMR in the vicinity of the impurity. From an experimental point of view, this is 
achievable by using conductive Atomic Force Microscopy approach as realized for instance in 
the following reference [? ] where the authors mapped out the spatial variations of the I(V) 
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characteristics through a tunnel barrier. From a theroretical point of view, a theory of local 
impurity assisted tunnelling in MTJ was recently developed jsl. Tight binding model and 
Kubo formalism were used to calculate the spin-dependent tunnel current through the MTJ. 
In this earlier paper, the I-V characteristics were not investigated in detail. Furthermore, 
the dependence of spin-dependent current on the position of cross section plane relative to 
the position of impurity was not calculated. 

In the present paper, we report on a theoretical study of the spatial distribution of 
spin-dependent current across the plane of a magnetic tunnel junction. The local I-V char- 
actertistics as well as the local TMR amplitude are calculated for a single impurity and for 
a random planar distribution of impurities inside the barrier. In this theory. We adopted 
the free electron model with exchange splitting for the ferromagnetic electrodes and used 
the nonequilibrium Keldysh technique jOf to calculate the transport properties which are 
nonlinear functions of the applied voltage. 

The MTJ is described as a three layers system, consisting of two thick ferromagnetic 
electrodes F separated by an insulating layer, I. Inside the barrier, a single nonmagnetic 
impurity with attracting potential is located at a given distance from the F/I interface. 
The two cases of parallel and antiparallel orientations of the F-layers magnetization were 
investigated. 

The F-electrodes are connected to the reservoirs with chemical potentials /ii and /i2 so 
that /i2 — /^i = eV, where V is the applied voltage. 

To calculate the current through the system, the Keldysh Green function G^^ and ad- 
vanced and retarded Green functions G"^ and G^ must be calculated. By solving the Dyson 
equation, we found that 



Go^(r,ro)WG'o-+(ro,r') , Go+lr, ro)H^G^(ro, r') 

(1) 



l-l^Go«(ro,ro) 1 - iyG^(ro, ro) 

, Go«(r,ro)W^G'o-+(ro,ro)iyG^(ro,r') 



(1 - W^G«(ro, ro)) (1 - WG^{v^, ro)) 

where Gq ^(r, r'), ^^(r, r') and Gq^TjT') are the Green's functions for the system in the 
absence of the impurity and the potential of the impurity V was represented as a 5-function: 
y(r) = Wal6{z — zq)6{p-Po), tq = {po, zq) is the position of the impurity, is its effective 
radius, W is its amplitude. The explicit expressions for G^, G^, G ^ have the following 
forms: 
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G?(r, r') = / d'K^-^^^^ {E{z„ z) [q{z,) + rk,] + E-\z,, z) [q{z,) - rk,]} 
X {E{z\ z^) [q{z^) + ik^] + E-\z', z,) [q{z^) - ik^]} , 

G^{r, r') = / 'i^'^fvTTn^ ^^^'^^ [^("^^ " '^'^ + ^"'(^2' ^^^^^^ + '^'^^ . ^ 

2^ q{z)q{z')den* /gN 

X 2i) - iki] + ^i) [g(^i) + iki]} , 

Go-+(r,r') = Id'^'^^!^^f^^^^{Eiz',z,) [qiz,) + ^k,] + E-\z' , z,) [q{z,) - ^k,]} 
X Z2) [qiz2) - ik2] + E~\z, Z2) \q{z2) + ik2]} 

(4) 

X - z/ci] + E-\zi,z) [q{zi) + i/ci]} , 

where 



fcl = ^/i^(£-Al)-'^^ 



A;2 = ^i^(e-A2 + eV)-«:2, 

den = {E{zi, Z2) [q{z2) - ^^2] [q{zi) - iki] - E^'^^zi, Z2) [q{z2) + ik2] [q{zi) + iki]} , 

K, is the electron momentum perpendicular to the plane of structure, e is the energy, zi and 
Z2 are the positions of the F/I interfaces, Ai and Ai denote the positions of the bottom of 
the energy band for spin up and down subbands. 

riL = /^(e) and = f^{e + eV) are Fermi distribution functions in the left and right 

^22 

reservoirs and height of potential barrier above Fermi level. 

In (III),©,® and (jH) p and z are in the plane and perpendicular to the plane coordinates, 
and we consider that z and zq are situated within the barrier. We have to take into account 

t I t I 

that all Green functions are matrices in spin space. k{p, k^p, kJ^p, k2p are Fermi wave vectors 



4 



4.0- 



3.5- 
3.0- 
2.5- 

7^ 

2.0- 

N 
II 

N 

3 1.5- 
1.0- 
0.5- 
0.0- 



— ■— Configuration tt Spin up 

— " — Configuration ti Spin up 

— » — Configuration ti Spin down 

— •—Configuration tt Spin down 



^„=5 (A) 



P(A) 



FIG. 1: Dependence of the current for different spin channels and P and AP configuration on the 
distance from the impurity in the plane of the structure at z=15 A. k^p = 1.1 A ^ , kp = 0.6 A ^, 
qo = 1.0 A~' 



of electron with spint (|) in the left and right F-electrodes. The current was calculated, 
using the following expression: 

eh f f dG'+{z,p;z',p) dG-+{z,p;z',p)\ 
MP-^)-l^J"e[ 5p 5j j^^^^ (5) 

In Fig^ and FigEl the dependencies of the currents in different channels (up and down 
spin) on coordinate p — Pq a.t one interface I/F {z2 = ISA) (another interface is at zi = 0) 
and inside the barrier at 2 = 10 are shown. In this calculation, the impurity is assumed to 
be positioned at po = and zq = SA. 

In the vicinity of the impurity, a hot spot of radius approximately equal to 6A may 
be observed. The current density in the center of the hot spot exceeds the value of the 
background current by several orders of magnitude. In FigOl the TMR dependence on the 
distance from the impurity at different z is shown. It is interesting that the value of TMR 
in the vicinity of the impurity exceeds its background value (TMR for the ideal structure is 
equal 0.013) by more than an order of magnitude. Furthermore, in some cases, regions of 
p — Po exist in which the TMR becomes negative. 

FigH shows the I-V characteristics for positive and negative applied voltage. These curves 
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FIG. 2: The same dependence at z = 10 A. 
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FIG. 3: Dependence of TMR on the distance from the impurity in the plane of the structure at 
different z. For parameters see FigQ 



are quite asymmetric with respect to the sign of the voltage. This asymmetry is related to the 
asymmetric positioning of the impurity inside of the barrier. It is particularly pronounced 
if the potential of the impurity is chosen so that the bound (resonance) state of electrons 
with spin up is located near the Fermi energy for the positive applied voltage = 1.2 V, and 
if this bound state lies below the Fermi energy for negative voltage. This diode behavior 
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FIG. 4: Local I-V curve at p = po and z = 15 A for the case of single impurity. 
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FIG. 5: I-V curve in the case of the layer of impurities at zq = SA and x = 0.5. 

was demonstrated so far in the case of a single impurity. We next investigate the case of a 
finite concentration of impurities. 

In this case, we consider the same magnetic tunnel barrier structure with a monolayer 
of impurities of finite atomic concentration x, situated closer to one of the F/I interfaces. 
To solve the problem, as a first step, we have to find the coherent potential and effective 
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Keldysh Green function G^^^. By solving the Dyson equation in the Keldysh space, the 

'TT 

Go~^ {z, zo)T,'^Gq {zo, z') Gq {z, zq) J]^Gq~^ {zq, z' 



following expression was obtained for G^ 



G-^{z,z') = G,+ iz,z') + 



Gq{z,zo)I1 ^Gq{zo,z') Gq {z, zo)T.^Gq^ (zq, zo)T.^Gq (zq, z' 



(1 - G^ (^0, ^o) S^) (1 - (^0, ^o) S^) ^ (1 - (^0, ^o) S^) (1 - (^o, ^o) S^) ^ ^ 

where T,^^'^^ are the coherent potentials (CP.) for the retarded and advanced Green func- 
tions, which have to be found from the C.P.A equation: 

t = (i-x) ~ + (x) ^^"^ ~ = (7) 

1 - (e^ - ^)Gesizo, po; zq, po) 1 - (e^ - T.)Gcs{zo, po; zq, po) 

where and are the onsite energies of the host {AI2O3) and the impurity (Al) and 
S-+ = |K + n^)(S^-S^). 

Now to calculate the I-V characteristic, we can use the previously found expression for 
Gaa^^^^^ and substitute it into the expression 

In Figl^l the I-V characteristic in the AP configuration is shown. An asymmetry of the 
curve on the sign of the applied voltage is clearly visible. 

Such a structure may be prepared for instance by sputtering a thin layer of Al on the 
bottom F-electrode, then oxidise it in Alumina. Thenafter, a second thicker layer of Al 
is sputtered on the already formed Alumina barrier but this second layer is subsequently 
underoxidized so that a thin layer of the random alloy AlxAl203(i-x) remains inside the more 
or less ideal insulator AI2O3 at an assymmetric location within the barrier. 
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